Unscreened Coulomb repulsion in the one dimensional electron gas 
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A tight binding model of electrons interacting via bare Coulomb repulsion is numerically investigated 
by use of the Density Matrix Renormalization Group method which we prove applicable also to very 
long range potentials. From the analysis of the elementary excitations, of the spin and charge 
correlation functions and of the momentum distribution, a picture consistent with the formation of 
a one dimensional "Wigner crystal" emerges, in quantitative agreement with a previous bosonization 
study. At finite doping, Umklapp scattering is shown to be ineffective in the presence of long range 
forces. 



One dimensional electron models are often used 
to interpret the behavior of strongly anisotropic 
physical systems in condensed .-.matter, like organic 
conductorsa, charge transfer saltsa and certain semicon- 
ductor nanostructuresu. These systems can be modeled 
in terms of a tight binding hamiltonian with effective elec- 
tron repulsion, which can be either short or long ranged, 
depending on the extent of interchain screening. The 
prototype of the short range models is the exactly solv- 
able one dimensional repulsive Hubbard modelcl which 
is known to be in the universality class of Luttinger 
liquidstl, showing metallic properties and antiferromag- 
netic spin correlations at every finite doping. However, 
this picture may dramatically change if screening is not 
effective in reducing the range of the bare Coulomb po- 
tential. This possibility is well known in quantum chem- 
istry where the Pariser-Parr-Pople model (PPP), describ- 
ing conjugated polyenes^, has exactly the same struc- 
ture of a tight binding model with unscreened Coulomb- 
like interaction. This hamiltonian has been studied in 
small lattices_by variational methods, like unconstrained 
Hartree-FockQ, suggesting the development of a charge 
density wave (CDW) leading to an antiferromagnetic in- 
sulating ground state in the thermodynamic limit. Exact 
diagonaliaations, recently performed in systems with few 
electronsa, instead provide evidence in favor of a metallic 
behavior. Unfortunately, finite size effects inhibited the 
study of correlation functions leaving open the two pos- 
sibilities of a Luttinger liquid or of a CDW metal where 
the charge carriers can be identified as sliding density 
waves, while electron-like quasiparticles are absent. The 
latter scenario was in fact proposed by Schulza in a sem- 
inal bosonization study. Here, the picture of a one di- 
mensional Wigner crystal emerged, characterized by ex- 
tremely long range tails for the charge correlation func- 
tions at the density dependent wavevector 4&p, together 
with weaker antiferromagnetic correlations at 2kp. 

In order to clarify the physics of one dimensional sys- 
tems with unscreened Coulomb repulsion, we performed 
a Density Matrix Renormalization Group (DMRG) study 
for the PPP model. Results for the excitation spectrum 
in the spin and charge channels are presented. The den- 



sity and magnetic structure factors have been computed 
both at half filling and in the doped system. A physical 
insight on the nature of charge carriers has been obtained 
through the study of the momentum distribution of the 
electrons. 

DMRG is an extremely accurate numerical method, es- 
pecially devised for one dimensional problems, which can 
easily handle system sizes two or three times larger than 
usual diagonalization algorithms, thereby drastically re- 
ducing finite size effectsta. DMRG gives the exact spec- 
trum of the hamiltonian in a reduced Hilbert space and 
therefore it provides a variational bound to the exact 
ground state energy. The method also allows for a self 
consistent evaluation of the errors introduced by trunca- 
tion, by checking the unitarity sum rule satisfied by the 
exact density matrix, which reflects the completeness of 
the Hilbert space. 

The model we have studied is defined by N electrons 
on a L site ring with hamiltonian 
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is a fermionic annihilation operator, rii — 
J2a c la c i,<? i s the density operator at site i with aver- 
age value n = N/L and d%j is the chord distance on the 
circle: 
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The parameter 70, which controls the strength of the 
long range Coulomb repulsion, is fixed at 70 p,l. 053907, 
following the Mataga-Nishimoto prescriptionlllL Here we 
present extensive results for a representative value inter- 
action strength U/t = 13.55. 

The physics underlying the half filled case (i.e. n = 1) 
is well established. In the strong coupling limit (i.e. 
U — ► 00), charges are frozen and the low energy con- 
figurations just correspond to different spin orientations: 
in a bipartite lattice, the ground state is always a non 
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degenerate spin singlet and it is expected to show cor- 
relations typical of one dimensional antiferromagnets: 
< S n ■ So >^ (— l) n /n, like in the Heisenberg model. 
Spin excitations are gapless with linear dispersion re- 
lation. Analogously to the Hubbard model, for every 
strength of the Coulomb repulsion, a charge gap develops, 
the system flows toward strong coupling and its correla- 
tion functions have the same structure as in the U — > oo 
limit. The open questions regard the effects of doping 
in the presence of long range repulsive interactions. In 
the Hubbard model, as soon as we insert holes (or extra 
electrons) into the lattice, the charge gap closes and both 
the charge and the spin branch of the excitation spectrum 
acquire linear dispersion relation, characterized by finite 
charge and spin velocity. In this class of spin isotropic 
Luttinger liquids, the long wavelength properties of the 
system are determined by a single dimensionless param- 
eter K p governing the power law decay of all correla- 
tion functions. The antiferromagnetic correlations in the 
ground state are characterized by the wavevector q = 2fcp 
and decay as x~^ 1+Kp ' while the density correlations are 
dominated by oscillations at q = Akp (i.e. wavelength 
A = 1 jn) and behave as x~ AKp . The momentum dis- 
tribution has no longer a sharp jump at the electronic 
Fermi momentum but maintains a power law singularity 
of the form ±|q-fc F | Q with a = (K p + 1/Kp - 2)/4. The 
analysis based on bosonization techniquesO showed that 
the long range nature of Coulomb interaction is expected 
to change this picture: while the spin sector is unaffected 
by interactions, the charge spectrum now becomes non 
analytic: uj q ~ qy/\ hi <?|- Formally this corresponds to 
the K p — > limiting case. In fact, the leading power law 
behavior of spin correlations is just cos(2kpx)/x while 
charge correlations decay more slowly than any power 
law and the momentum distribution is continuous with 
all its derivatives at the Fermi momentum. Actually, the 
bosonization study allowed to determine the precise form 
of the asymptotic behavior, which turns out to be: 
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with c and c' non universal positive constants. However, 
Umklapp scattering, possibly leading to CDW instabil- 
ities, has been neglected in the bosonization analysis, 
as remarked in RefH. Other possible scenarios involve 
charge ordering with different periodicity, i.e. the one di- 
mensional analog of a "stripe" phase, or the formation of 
a gap in the spin excitation spectrum, as a precursor to 
phase separation, which is inhibited by long range forces. 

In order to understand which picture correctly de- 
scribes the physics of correlated electron systems with 
long range interactions, we numerically studied the PPP 
model (|l]) at the two electron densities n — 1/2 and 
n = 3/4 corresponding to kp = tt/A and kp = 3n/8 re- 



spectively. These choices of filling factors are suggested 
by the necessity to perform accurate size scaling keeping, 
at the same time, a limited number of degrees of free- 
dom. Calculations have been performed in lattices up to 
L = 80 sites and total number of electrons up to N = 60. 
The chosen electron densities are represented by small 
fractions, hence we expect that commensurability effects 
may be enhanced for these cases: If the model is prone 
to a CDW instability, the systems studied in this work 
should clearly suggest the tendency toward charge order- 
ing. Our DMRG code is written in such a way to cope 
with long range potentials. The dimension of the reduced 
Hilbert space is always larger than 10 6 , the truncation 
error is at most 10 -5 and the correlation functions are 
translationally invariant up to a relative error of 2 x 10p^. 
Further details on the algorithm can be found in Ref.lia. 

As a first step, we calculated the charge spectrum, i.e. 
the ground state energies E(N + 1) and E(N — 1) ob- 
tained by adding and removing one electron to the ref- 
erence state with N — nL particles. In order to have 
a non-degenerate singlet ground state, we imposed peri- 
odic or anti-periodic boundary conditions thereby realiz- 
ing the closed shell condition in the non-interacting limit 
on the reference state. The finite size gap is defined as 
the difference between the upper and lower estimate of 
the chemical potential: 



A, 



- [E(N + 1) + E(N - 1) - 2E(N)] 
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Although A p does not coincide with the true charge gap 
we can infer that the charge gap is zero if A p vanishes 
in the thermodynamic limit. Conversely, the spin gap 
A a- is just the energy difference between the singlet and 
the triplet spin sectors at fixed number of particles. The 
finite size scaling of our results, shown in Fig. 1, clearly 
indicates the gapless nature of the excitation spectrum, 
ruling out the possibility of charge ordering in the ground 
state and confirming the irrelevance of Coulomb repul- 
sion on the structure of spin excitations in agreement 
with previous diagonalization dataa. The analysis of cor- 
relation functions provides a deeper information on the 
physical nature of the ground state. In particular, we 
studied the charge and spin structure factors defined by: 

R(q) = = ~ E eiqil ~ m) < ( n < ~ ") " ") > 

l,m 

S(q) = = ±J2 em ~ m) < S ? SZ m> (5) 
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The Fourier transform R(q) of the equal time density- 
density correlation function shown in Fig. 2 displays a re- 
markable collapse of data relative to different sizes, with 
the single exception of the peak region at q = Akp. The 
peak value, plotted as a function of L in Fig. 3, is well fit- 
ted by the expected form R(ikf) ~ L exp(— 4cVhi L) de- 
duced from the bosonization analysis (0) suggesting that 
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quasi long range order develops in the system. Notice the 
appreciable curvature of the data which reflects the pres- 
ence of the exponential term in the fitting formula and, in 
turn, indicates a very slow decay of density correlations 
in real space. The holes effectively repel each other and 
stay almost rigidly at the maximum attainable average 
distance A = 1/n but, as expected, hole correlations van- 
ish at large distance, in agreement with Eq. (||). A much 
weaker singularity seems to be present at wavevector 2k p, 
but the data do not allow for a systematic analysis of 
this further feature, also predicted by the bosonization 
analysis of RefB. Other interesting information can be 
extracted from the small q behavior of R(q) which is re- 
lated to the low energy excitation spectrum. According 
to the bosonization formulas, the structure factor should 
behave as R(q) ~ q\ lng| -1 / 2 . This form well represents 
the DMRG data, as shown in Fig. 4. Note that such an 
expression is quite specific to the bare Coulomb interac- 
tion and differs from the usual Luttinger liquid results, 
where R(q) ~ K p (q/iv). A similar analysis can be car- 
ried out for the magnetic properties of the model. In 
Fig. 5 the magnetic structure factors S(q) correspond- 
ing to the previously discussed choices of parameters are 
shown. The collapse of points on the same smooth curve 
should be again appreciated. The dominant peak oc- 
curs at q = 2fcp, as expected, and its height, reported 
in Fig. 6, scales as ao — ai(v / lnL + 1/c) exp(— cVhiL) 
in agreement with the bosonization prediction (^|). Here 
ao and a\ are fitting parameters, while the constant c is 
fixed to the same value obtained from the fit of the den- 
sity structure factor shown in Fig. 3. Contrary to the 
results obtained for R(q), the small q behavior of S(q) 
is accurately given by the linear relation S(q) ~ (q/Air) 
corresponding to K a = 1, as expected for all the gapless, 
SU(2) invariant, one dimensional electron systems. 

The previously discussed properties of the PPP model, 
emerging from the DMRG analysis, are consistent with 
the Wigner crystal model suggested by the bosonization 
study: away from half filling, the system is gapless both 
in the charge and spin sectors but, at the same time, 
it supports extremely long ranged density correlations. 
Following Ref.l, we may consider two different pictures, 
consistent with the observed metallic properties of this 
model. The first one corresponds to a Luttinger liquid 
scenario, where the electrons have a sort of "Fermi sur- 
face" satisfying the Luttinger theorem. In this case, al- 
though the low energy excitations are always collective in 
one dimension, the charge carriers in the system may be 
still identified with dressed electrons. Instead, accord- 
ing to the second picture, the absence of a gap in the 
charge spectrum is due to the Wigner crystal nature of 
the ground state and to the long range decay of density 
correlations. The quasiparticles cannot be identified with 
electrons any more and consequently the electron mo- 
mentum distribution is smooth at the Fermi momentum, 
with an extremely weak essential singularity, as reported 
in Eq. (|3|). DMRG data for the momentum distribution 
for given spin n(q) are shown in Fig. 7, while a size scal- 



ing of the finite size jump at q = kp is plotted in Fig. 8. 
Again, the asymptotic scaling 

An(fc F )~Le~ c '( lnL > 3/2 (6) 

predicted on the basis of the bosonization result (||) gives 
a quite good representation of numerical data. 

In conclusion, we have tried to clarify the physics un- 
derlying the presence of unscreened Coulomb interaction 
in one dimension. Accurate DMRG data are fully consis- 
tent with the picture of a Wigner crystal showing metal- 
lic properties, with slowly decaying density correlation 
and gapless charge and spin excitation spectra. At least 
in one dimension, this numerical investigation rules out 
other scenarios, like pinned CDW or even some precursor 
of off-diagonal long range order. 
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FIG. 1. Charge gap A p and spin gap A CT , in units of t, 
versus the inverse chain length L. The density is n = 3/4. 



1 -- 



cc 0.25 - 



o L = 


56 


□ L = 


48 


oL = 


44 


A L = 


40 


< L = 


36 


v L = 


32 


oL = 


28 


+ L = 


24 


»L = 


20 


* L = 


16 




FIG. 2. FIG. 2. Fourier transform R(q) of the equal time 
density-density correlation function. The density is n = 1/2 
in (a), n = 3/4 in (b). 



0.8 



16 



24 



32 



40 



48 



56 



FIG. 3. Fit of R(4kf) with respect to the bosonization 
result. The circles are DMRG results for density n = 1/2 and 
the solid line is a function of the form aL exp(— 4c\/ln L) with 
a = 1.93 and c = 0.54. 
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FIG. 4. Small momentum behavior of the charge structure 
factor. We plot R(q) divided by the expected limiting form 
q\ ln(q)|~ 1/2 . The density is n = 3/4. 




FIG. 5. Fourier transform S(q) of the equal time spin-spin 
correlation function. The density is n — 1/2 in (a) and 
n = 3/4 in (b). 




FIG. 7. Momentum distribution function for given spin 
n(q). The density is n — 1/2 in (a) and n — 3/4 in (b). 
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FIG. 6. Fit of S(2kj) with respect to the bosonization 
result. The circles are DMRG results for n — 1/2 and the solid 
line is a function of the form a a — a i ( VlnL+ 1 /c) cxp ( — cVln L) 
with oo = 1.12, ai — 0.68 and c = 0.54 (obtained from Fig. 
3). 
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FIG. 8. Fit of the finite size jump at kf of the momentum 



distribution with the expression: An(kf) = aLe 



-c' ln(L) 3 / 2 



The parameters are: a = 0.493781, c' = 0.631772 for n = 1/2 
(curve A) and a = 0.491581, c' = 0.611964 for n = 3/4 (curve 
B). 
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